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A general framework is developed to investigate the properties of useful choices of sta¬ 
tionary spacelike slicings of stationary spacetimes whose congruences of timelike orthog¬ 
onal trajectories are interpreted as the world lines of an associated family of observers, 
the kinematical properties of which in turn may be used to geometrically characterize 
the original slicings. On the other hand properties of the slicings themselves can directly 
characterize their utility motivated instead by other considerations like the initial value 
and evolution problems in the 3-plus-1 approach to general relativity. An attempt is 
made to categorize the various slicing conditions or “time gauges” used in the literature 
for the most familiar stationary spacetimes: black holes and their fiat spacetime limit. 
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2 D. Bini, E. Bittencourt, A. Geralico and R.T. Jantzen 

1. Introduction 

A century after the birth of general relativity, we now take for granted the existence 
of various stationary spacelike slicings of stationary spacetimes which have certain 
special geometrical properties useful in studying the astrophysical consequences of 
say, black hole spacetimes. Many of these slicings arise from geometrical properties 
of their irrotational congruences of orthogonal timelike trajectories, interpreted as 
the world lines of an associated family of observers which may be either geodesic 
or accelerated. Other slicings are instead characterized by the intrinsic or extrinsic 
geometry of the slicing itself. We here survey the various categories of such useful 
slicings for nonrotating and rotating black hole spacetimes, but starting with the 
limiting flat Minkowski spacetime which allows the greatest variety of examples of 
special slicings. 

Among the observer-defined slicings of black hole spacetimes is the Boyer- 
Lindquist time coordinate slicing associated with the usual stationary accelerated 
observers referred to equivalently as fiducial or locally non-rotating or zero angular 
momentum observers, abbreviated as FIDOs, LNOs or ZAMOs. Geodesic observer 
families instead characterize the rain, drip and hail coordinate systems m which 
include the Painleve-Gullstrand coordinates [314] and their generalizations For 
nonrotating black holes the latter are also characterized by the intrinsic curvature 
properties of their associated slicing, whose induced geometry is flat. Smarr and 
York pioneered linking the preservation of kinematical properties of the slicing 
to the choice of time lines in an evolving spacetime, while constant mean curvature 
slicings were seen as privileged from the point of view of the initial value prob¬ 
lem even earlier [HHE]. Special spherically symmetric slicings of the nonrotating 
Schwarzschild black hole spacetime were first considered by Estabrook et al nnm]. 
Recent work has investigated the geometry associated with analogue black holes 
PTT^ and shown that neither intrinsically flat nor conformally flat slicings of the 
Kerr spacetime exist |14l15l1til5) . 

On the other hand, the analysis of the Gauchy problem of general relativity has 
also led to the introduction of spacetime slicings useful in simplifying the evolution 
equations, like harmonic slicing [17] (and closely related time gauges [18]), defined 
by requiring the time coordinate of the spacetime associated with the slicing to be a 
harmonic function. Eurthermore, the numerical relativity study of multi-black hole 
dynamics m takes advantage of the use of “hyperbolic slicings,” requiring spatial 
compactification techniques at infinity [20121122] . as well as horizon penetrating 
coordinates like Painleve-Gullstrand coordinates. We refer to these as “analytic 
slicings,” belonging to the Gauchy problem literature, in contrast with those of a 
“geometrical” nature. 

Motivated primarily by a desire to better understand the underlying geometrical 
structure of these spacetimes, we systematically review, develop and discuss special 
slicings of black hole spacetimes together with their flat Minkowski limit. We use 
geometric units with c = 1 = G. Greek indices run from 0 to 3 and refer to spacetime 
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Slicing black hole spacetimes 3 

quantities, while Latin indices from 1 to 3 and refer to spatial quantities. 

2. Spacetime, general coordinates and spacelike 3-surfaces 

Let be a generic set of coordinates adapted to a slicing of spacetime by 

spacelike hypersurfaces of constant values of the time coordinate t and write the 
metric as 

= ga,ada:“da;^ = gttdf^ + 2giadtda;“ + gafcdx^da;^, (a, 6 = 1,2,3). (1) 

The convenient Wheeler lapse-shift notation |23| re-expresses the metric in the form 

d- gab{Ax‘^ + lV“dt)(dx^ -k N^dt) 

= -N'^df + , ( 2 ) 

defining the lapse function N and shift vector field which satisfy 

gtt = -{N^- N,N -), gta = Na = gabN^, (3) 

while the 1-forms {oJ°'} are orthogonal to the unit timelike 1-form = — 7Vd|3 
associated with the unit normal vector field n to the slicing 

n = j^{dt - N-da). (4) 

Here we use the equivalent shortened notations for the the spatial coordinate frame 
da = d^a = Ca- The induced Riemannian 3-metric on the time slices is simply 

^^^ds^ = gabdx^'dx '^, (5) 

with spacetime volume element ^/—g = N^det{gab) ■ Note that {N dt,uj°‘} is the 
dual frame to the frame {n,da} reflecting the orthogonal decomposition of each 
tangent space adapted to the slicing and its normal direction. 

For a stationary spacetime, one can choose these coordinates so that the station¬ 
ary symmetry corresponds to translation in the time coordinate t, with associated 
Killing vector field dt- To transform to another stationary slicing, without loss of 
generality one can consider restricted choices of the new spacelike time coordinate 
of the form 

T = t-f{x^,x^,x^) (6) 

which retain the time lines of the original coordinate system if the spatial coordi¬ 
nates are not changed. One is still free to choose new time lines by changing the 
spatial coordinates as well, but unless the time lines are associated with a Killing 
vector field, the metric will become explicitly time-dependent. 

On a generic slice E of the new slicing, the 1-form 

dT^dt-fadx\ fa = daf, (7) 


The symbol denotes here the completely covariant form of a tensor A. 
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4 D. Bini, E. Bittencourt, A. Geralico and R.T. Jantzen 
vanishes identically, i.e., 

dT|j, = 0 —>■ dt = /adx“ (on E). (8) 

If we retain the spatial coordinates and only introduce this new time coordinate, 

t = T + = (9) 

then one must distinguish the spatial coordinate frame vector fields tangent to the 
old and new time coordinate hypersurfaces 

dT=dt, Ea=dx<^ =da + fadt, dX“ = dx“ . (10) 

Re-expressing the spacetime metric then leads to 

(4)ds2 = + hab{dX^ + L^dT){dX^ + L^dT ), 

=-L'^dT^ + (11) 

where hab (the induced metric on E) and L and Lc (the new lapse and shift) are 
given by 

^ab — 9ab d~ (Jtafb d~ fjtbfa d~ Qttfafb : 

Lc = gtc + gttfc ) L‘^ = h‘^‘^Ld , = —gtt + . (12) 

Similarly to Eq. we then have 

QTT — (-^ LqL ) — 9tt-! 9Ta — Ld — hdbL — gia Qttfa ■ 

One can also evaluate the contravariant metric 

1 

TT ^ Ta ^ lab ab /-i a\ 

g =~J 2 ^ 9 =J 2 ^ h =g (14) 

while the spacetime and spatial metric determinants satisfy 

^/^=LVh, h = det{hab)- (15) 

Note that {LdT, 17“} is the dual frame to the frame {V, Eq} adapted to the 
orthogonal slicing decomposition of the tangent space and that on E one has | ^ = 
dX“ = da;“. Any “spatial tensor” has only Latin indexed components allowed to 
be nonzero in this frame. 

The unit timelike 1-form normal to the T = const slicing is given by 

Af'’ = -LdT = -L{dt - /adx“), (16) 

with associated unit timelike normal vector field 

Af=^{dT-L-Ea). (17) 

In turn Af can be expressed in terms of n as 


Af = 7 (A/’, n)[n + v{Af, n )], 


(18) 
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Slicing black hole spacetimes 5 

where is the relative velocity of J\f with respect to n and j{Af,n) the 

associated gamma factor, explicitly 

7(Ar, n) = ^{l + N<^h ), nf = , (19) 

as follows from Eq. (HH) after re-expressing dt and da:“ in terms on n*' and a;“. A 
straightforward calculation shows that the new lapse function and shift vector field 
for the same time coordinate lines are given by 

L = 7 (Ar, n)N[l - KAT, n),N<^] = , 

L = L^da = L 7 (Ar, n) n)“^ da • (20) 

The above decomposition gives a more transparent kinematical meaning to the 
various quantities, as we will show below. 

Starting from the spacetime unit volume 4-form rjfaaP'y, one can associate with 
any timelike unit vector field U, whether n or Af, a spatial volume 3-form ri(U)ap~i = 
U^rifaap-y which can be used to define the cross product xjj and the curl operator 
curlj/ in the local rest space LRSu of U, as well as a spatial duality operation for 
antisymmetric spatial tensor fields. In a spatial frame adapted to that subspace, for 
spatial vector fields X and Y in that subspace, one has 

[X XuY^ = r,{U)\aX’^Y\ [cmlu X]‘^ = r,iU)\a^{Uf X% (21) 

where the spatial covariant derivative V(C/) of any tensor field (including U) is 
obtained by projecting all indices of the spacetime covariant derivative of that 
tensor into the local rest space using the associated projection tensor P{U)°‘p = 
(5“/3 -f U°‘Up whose fully covariant form P{U)^ is the spatial metric. 

We conclude this section by introducing the relevant tensor quantities needed 
to evaluate both the intrinsic and extrinsic curvature of a typical slice S as well as 
provide a geometrical characterization of the kinematical properties of its normal 
congruence M. 

(1) Intrinsic curvature of S. 

This is obtained evaluating the Riemann tensor components R{h)abcd of the 
3-metric h induced on E, i.e., 

= hcddX^dX'^. ( 22 ) 

Note that in the three-dimensional case the Riemann tensor is completely de¬ 
termined by the associated Ricci tensor. 

(2) Conformal flatness of E. 

The Cotton tensor associated with the induced 3-metric (l22)l is given by 

[R{h)\b - lR{h)S‘^[b'^ 


C{hTba = 2 


;c , 


(23) 
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6 D. Bini, E. Bittencourt, A. Geralico and R.T. Jantzen 

where all operations including the covariant derivative f\a = V(A/')a/ refer to 
the 3-metric. A vanishing Cotton tensor characterizes the conformal flatness of 
the spatial metric. Taking the spatial dual of the two covariant indices gives 
the equivalent Cotton-York tensor 

= (24) 

which is symmetric because of the twice contracted Bianchi identities of the 
second kind 

C'(h)“^ = = Y{hY^ , (25) 

where 

Y{hr^ = = [ScurW i?(/i)]“' , (26) 

and the symmetric tensor curl “Scurl” operation [24125126] is defined by 

[ScurW R{h)Y^ = . (27) 

The spatial metric h is conformally flat if the spatial Ricci tensor R{h) has 
vanishing Scurl. 

(3) Extrinsic curvature of E. 

This is obtained evaluating the Lie derivative of the spacetime metric along 
the unit normal A/" to E and projecting the result orthogonally to J\f and raising 
an index to make a mixed tensor, with an extra numerical factor 

^(■ A ^)“/3 = - [£y 9]^s P{^)% • ( 28 ) 

Its trace Tr [K{J\f )] = K{MY p is the mean curvature of the slice. The constant 
mean curvature (CMC) time gauge is a slicing with constant mean curvature on 
each slice, though it may vary from slice to slice. A maximal slicing instead has 
vanishing mean curvature on every slice. When the tensor K(Af) itself vanishes, 
the slicing is called totally geodesic, or extrinsically flat. An invariant character¬ 
ization of the extrinsic curvature can be obtained by studying its eigenvalues, 
namely those of the 3x3 matrix of components {K{J\f)°‘b) in an adapted frame 
since it is a spatial tensor. The three eigenvalues in turn can be expressed in 
terms of the three scalar trace invariants of the powers of the extrinsic curvature 
Ti = Tr [K{N)\, Ti = Tr [K{MY\ and 7^ = Tr [K{M)% 

(4) Kinematical fields associated with M'. acceleration, vorticity, expansion and 
shear. 

These are obtained by decomposing the covariant derivative of J\f into its 
irreducible parts under a change of frame 

= -a{MYMp - uj{MYp + e{MYp , (29) 

namely the acceleration a(A/')“ = VjgAf‘^ of the congruence, its vanishing vor¬ 
ticity tensor uj{Af) = 0 since the congruence is hypersurface-forming and the ex¬ 
pansion tensor 0{J\f) = —K{M). The scalar expansion Tr [0(A/’)] = —Tr [K{M)\ 
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Slicing black hole spacetimes 7 


of the congruence is just the sign-reversed mean curvature of the slicing, while 
the tracefree part of the expansion tensor is the shear tensor. 


Finally, there exist slicings motivated more by analytic considerations than ge¬ 
ometrical ones. One example is represented by the so-called harmonic slicings. In 
this case the T = const foliation is characterized by a harmonic condition 


□T = = 0 . 


(30) 


A similar condition imposed on all spacetime coordinates specifies the so-called de 
Bonder gauge choice of coordinates. In the literature, there are also variations of 
this condition (“1 -|- log” slicing, etc.), which will not be explored here (see, e.g., 
Ref. [U]). Eq. (1501) is equivalent to 


1 dT 



dT \ 


= 0 , 


which becomes in a coordinate system with time coordinate T = 


(31) 


dg{V^9^^)=0. 


Recalling Eq. 0 the previous equation can be also written as 



(32) 


(33) 


which in turn can be expressed as an evolution equation for the new lapse function 
L. 


3. Minkowski spacetime 

As a simple example of the problem of finding special slices in a given spacetime, 
let us consider the flat Minkowski spacetime geometry in an inertial time slicing, 
with its line element written in spherical coordinates to compare later with a black 
hole spacetime in Boyer-Lindquist coordinates 

ds^ = -df + dr^ -f r^de^ -H sin^ 6 . (1) 


The t = const hypersurfaces are both intrinsically and extrinsically flat. Because of 
their simplicity, the Minkowski slicing examples are useful in better understanding 
the corresponding general relativistic situations considered below. 

Consider a new spherical symmetric slicing by a time function T = t + f{r), 
retaining the spatial coordinates (and hence the orthogonal time lines). The induced 
metric on a typical slice E is 

(hab) = diag (1 - /^, r^, sin^ O) , (2) 


with the new lapse and shift functions 
I 


L = 


x/T^’ 


La = -frS^a , 


1 -/) 




L, 


L 


( 3 ) 
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One can define the new spatial frame 


El — dr + fr dt, E 2 — dg , E 3 — d^ , 


(4) 


with dual frame 


^ dr, = de, = del). 


(5) 


In this frame the extrinsic curvature tensor has components 


with trace 


( 6 ) 

(7) 


The intrinsic curvature is characterized in three dimensions equivalently either by 
the Riemann or Ricci tensor, with nonvanishing components respectively 




'^frfrr E{h')r(l)r<p 

^ — fr sin^ 9 


E{Ji) Q (ftO (f) — 


sin^ Off 
1-/2 


( 8 ) 


and 

R{hYr 


2/r/rr 

7(1 - /2) 


R{h)\ 


fr 

r2(l - /2)2 


[rfrr + fr{l-ff)] 


R{h)\, (9) 


and the curvature scalar is 

+ /-(I - fr)] • (10) 

Finally, the Cotton-York tensor (l24l) is identically zero for this slicing, independent 
of the choice of /(r), ensuring the conformal flatness of any spherical slicing. 

The new time coordinate hypersurfaces T = const have unit normal 


=—L{dt — vdr), J\f = ^[dt + v dr), (11) 


with radial relative velocity v = fr (which must satisfy \v\ < 1) and associated 
Lorentz factor 7 = 1 /\/l — i/2. The associated observers moving orthogonal to 
the slicing follow radial trajectories which are ingoing for ly < 0 and outgoing for 
> 0 in comparison with the usual static observers {v = 0) following the original 
time lines. The simple transformation / -A —f flips the radial motion of the new 
observers. 

Let us consider some explicit examples. In order to deal with dimensionless 
quantities, we introduce a positive scaling constant with the dimensions of a length, 
say A. For simplicity the solutions of the various conditions on the slicing function 
/(r) below will be given modulo an overall sign subject to the initial condition 
/(O) = 0 for the sake of easy comparison. 





















September 17, 2015 0:19 WSPC/INSTRUCTION FILE kextr'final'grqc 
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(1) CMC: 

Let Tr[iL(A/")] = kjA = const, with k a dimensionless constant. Eq. (O 
then gives 

fr kr CA^ 


\/r^ 


= X, 


( 12 ) 


where C is a dimensionless integration constant, implying that 


hrr = {l + X^y 


— .27 h^rp. 


(13) 


The function /(r) is then given by 

/(r) = f . 

Jo 


(14) 


(2) Maximal: 

These slices are characterized by the vanishing of the trace of the extrinsic 
curvature tensor Eq. © which leads to 

fr = y = , , (15) 


Y^l + ’ 


modulo a constant chosen to be 1 with no loss of generality. This equation can 
be integrated in terms of elliptic functions, i.e., 

/(.) = -1 {f (^„cco» . T) _ 2K (T^ I , (16) 

where F{x, k) and K{k) are the incomplete and the complete elliptic integrals 
of first kind, respectively (such that F{TT,k) = 2K{k)). Note that in the limit 
of r —>• oo we have /(r) —>• AK{l/y/2) r; 1.854 A. 

(3) Vanishing Ricci scalar. 

Requiring that the spatial Ricci scalar Eq. m vanish leads to a parabola 
of revolution about an inertial time axis orthogonal to its symmetry axis 

/(r) = 2yi(]4+7) - 2A, r=—^=^=. (17) 

v'l + rjA 

The induced metric is 

(/j-afi) = diag ^ ^^^ ,r^,r^sin^6l^ , (18) 

with extrinsic curvature 

(K(AA)%) = yjdiagQ,-l,-l) , V^[K{N)] = -\^, (19) 

and intrinsic curvature specified by 

(i?(/i)“6) = ^diag . 


( 20 ) 
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(4) Hyperboloidal: 

Smarr and York j6] introduced a special constant mean curvature slicing of 
Minkowski spacetime by translating a single sheet of a spacelike pseudosphere 
(hyperboloid) of a given fixed radius (and therefore fixed intrinsic and extrinsic 
curvature) along the time lines of an inertial coordinate system, in contrast 
with the geodesically parallel family of all such pseudospheres of varying radii 
centered on one spacetime point. In the metric ([IJ in spherical coordinates in 
an inertial coordinate system, consider the time slices T = const where T is 
given implicitly by 


{t + A-Tf-r^= 


( 21 ) 


so that /(r) = \/r^ + A^ — A, where A is the radius of the pseudosphere and 
we have chosen only one of the two possible signs, i.e., the plus sign. One hnds 
then the corresponding radial velocity given by 


n = 



( 22 ) 


The normal congruence to the T = const slicing defines the associated Smarr- 
York observers. The intrinsic metric of the slices evaluates to 

(hab) = diag ’ (^3) 

while the extrinsic curvature tensor 

K{Af)ab = -\hab, Tr [K{Af)] = - |, (24) 

and Ricci curvature tensor 

R{h)ab = ^Kh ■, Rih) =, (25) 

are both (spacetime) constant multiples of the spatial metric reflecting the 
constant intrinsic and extrinsic curvature conditions. 

The lapse function and shift vector field are given by 

/>^2 _ rv^ I 

, L = L°-da =I + ^ dr , (26) 

respectively. The spatial metric then takes the form 

rl-r2 

(3)ds2 = ^+r^dn^, dfl^ = d6»2 + Sin^ 6 d4i^ . (27) 

or transforming the radial coordinate by r = Asinhr]^ 

(3)ds2 = A'^idr)'^ + sinh^T? dfl^), (28) 



more familiar from cosmology. Note that when A —> oo the slice tends to 
be both intrinsically and extrinsically flat. This shift vector field satisfies the 
minimal-distortion equation of Smarr and York. 


The behavior of the slicing function /(r) as well as of the associated spatial 
velocity i/ as a function of r/A is shown in Fig. [T] in all cases discussed above. 
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Slicing black hole spacetimes 11 




Fig. 1. Minkowski spacetime. The radial behavior of the slicing function f{r) (rescaled by A) 
and the associated spatial velocity u = fr of observers normal to the slicing itself relative to the 
observers following the coordinate time lines are shown as functions of r/A in all cases discussed 
in the text: CMC, maximal, vanishing Ricci scalar and hyperboloidal slicings. The curve corre¬ 
sponding to the CMC slicing is drawn here for fc = 1; as fc approaches 0, this curve approaches the 
maximal slicing curve more and more. Note that although all these slicings have observer world 
lines which are initially radially outgoing at the origin of spatial coordinates, the CMC observers 
reverse direction to become radually ingoing. Only the hyperboloidal slicing is regular at the 
origin, while the remaining hypersurfaces have an asymptotically null conical singularity there. 
Reversing the sign of f{r) simply interchanges ingoing and outgoing directions for the relative 
motion of the new observers with respect to the inertial observers of the rest system. 


4. Schwarzschild spacetime 

Consider now the Schwarzschild spacetime representing a nonrotating black hole, 
whose line element written in standard coordinates (t, r, d, </)) is given by 

ds2 = - (^1 - dt^ +(^1- ' dr^ + r2(d02 + sin^ 0d<^2) _ 

Introduce an orthonormal frame adapted to the static observers following the time 
lines, i.e.. 


n = N-% , 

01 

II 

^0 : 
r 


(2) 

^ r sin 6/ 

with dual frame 

1 

II 

= N-^dr , 

A = rdd , 

= r sin 0 d4> , 

(3) 


where N = \/l — 2M/r. 

The t = const hypersurfaces form a slicing which is extrinsically flat (as the 
orthogonal hypersurfaces to the static Killing vector congruence dt), but not intrin¬ 
sically flat. The induced metric on these hypersufaces 

(3)ds2 = A - ' dr2 + r2(d02 + sin^ (4) 
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has nonzero Ricci curvature 

(i?(/i)“6) = ^diag (-2,1,1) , (5) 

with vanishing Ricci scalar R{h) = 0. 

Let us look for general slicings of the Schwarzschild geometry which are compat¬ 
ible with the Killing symmetries of the spacetime, i.e., spherical slices T = t + f(r). 
Their timelike unit normal is 

= — L{dt — fr dr), Af = ^{n + ef), (6) 

with relative velocity (see Eq. ®) 

V = V^ = N^fr ^ fr= N-^V (7) 

and associated Lorentz factor 7 = l/\/l — ly^. As in the Minkowski spacetime case, 
reversing the sign of f(r) interchanges the ingoing and outgoing radially moving 
observers associated with the new slicing. 

The induced metric on E is 


(hab) = diag ((yAf) 
and the new lapse function and shift vector field are 
L = Nj, La = -iyS\, = 

A (nonorthonormal) basis on E is given by 

El = dr + frdt ) E2 = de , E3 = d^ , 

with dual frame 


[i/ + w’’], = 

The acceleration a(Af) is given by 


1 




1 


= d 6 », = 


1 


y/ 


= d(j). 


a{N) = —dr{N^)n\ 

and the extrinsic curvature of the T = const slices then turns out to be 
K{Af) = - rX 0 -f sin^ 0 0 , 

where X = N^v, and the traces of its powers then have the following values 


Tr[iL(V)”] = (-1)" 


{drX)- + -A" 


(n = l,2,3). 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


The intrinsic curvature is described by the following nonvanishing components 
of the spatial Riemann tensor 

R{h)rere = = -^driN^) , R{h)e^ect = r^sm^0{l-^^N^), (15) 

Sin 0 iV 7 
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Slicing black hole spacetimes 13 

or equivalently by the spatial Ricci tensor, whose nonzero components are 

R{hYr = --N-fdr{N-f), R{h)%=RihY^ = lR{h)\ + ^{l--f^NY. (16) 

r Zr 

Finally, the spatial Ricci scalar is 

, (17) 

showing that r(l —= const corresponds to vanishing spatial scalar curvature, 
but jN = 1 to vanishing spatial curvature. 

The Cotton-York tensor (l24)l vanishes identically, so that spherical slicings are 
automatically conformally flat. 

Geodesic slicings correspond to a spatially constant lapse function L = Nj = 
Lq, as from Eq. (IT^ . In this case 

and the only surviving component of the spatial Riemann tensor is 

Rih)e4,0,j, = r'^ sitY 9 {1 - LI) , (19) 

whereas the spatial Ricci tensor is fully specified by 

R{hYe = R{hY^ = ^{l-Ll), ( 20 ) 

with R{h) = 2(1 — L'Y)lr^. The value of the parameter Lq determines the sign of 

the intrinsic curvature of the slices, corresponding to either vanishing (Lq = 1), 
positive {Lq < 1) or negative {Lq > 1) curvature. Note that interpreting Af as the 
4-velocity field of a family of (geodesic) test particles, the parameter Lq coincides 
with the (conserved. Killing) energy per unit mass of the particles 

E = -J\rt = Nj = Lo. ( 21 ) 

Therefore, Lq = 1 also represents the case of particles with vanishing radial velocity 
at spatial infinity; Lq > 1 can be associated with particles starting at infinity with 
nonzero (inward) velocity; finally, Lq < 1 corresponds to particles starting moving 
at a finite radial position. In the literature, adapted coordinates to these situations 
(in the special case of a Schwarzschild spacetime) are referred to as “hail” {Lq >1), 
“rain” (Lq = 1) and “drip” {Lq < 1) coordinates (see, e.g.. Ref. [2] and references 
therein). 

Let us consider some explicit examples. We will let C denote in each case the 
dimensionless integration constant appearing in the solution of the ordinary dif¬ 
ferential equation for the radial velocity all quantities being rescaled by the 
characteristic length scale of the background curvature (i.e., the mass M). For sim¬ 
plicity the solutions for the slicing function f{r) will be given up to an overall sign 
and the constant of integration will be chosen so that f{2M) = 0 if this function 
is finite at r = 2M, or lim^-s-oo /(?') = 0 if it approaches a finite limit at r —^ oo, 
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in order to compare the time slices from a common reference point at the horizon 
or at spatial infinity relative to the usual slicing. Typical behaviors of the slicing 
functions as well as of the corresponding linear velocities are shown in Fig. [2j 

(1) CMC: 

Let M Tr [K{J\f)] = k = const. Then using Eq. (imi with n = 0 one finds 

kr CM^ 


X = -- 


( 22 ) 


formally equivalent to Eq. (HU in the case of a flat spacetime, where C is a 
dimensionless integration constant, implying that 




and 


Kr = ^ = + X^y 


The function /(r) is then given by 

f{r) = f 

J21} 


X 


: dr . 


(23) 


(24) 


(25) 


I2M N^KN^+X^ 

Let k > 0. Then the radial velocity approaches jz —)> — 1 for large r, independent 
of the value of C. In contrast the behavior at the horizon depends on the sign of 
the quantity —8k+3C. In fact, in the limit r -A 2M one has n -A sgn(—Sfc+SC), 
i.e., 1 / —>■ 1 if C > |/c, while —1 if C < |fc. 

(2) Maximal: 

For k = 0 the above relations simplify to 

CM^ 


X = 


so that v = CM^KKr/K and 




2M 

r K 


and 




(26) 


(27) 


(28) 


which can be expressed in terms of elliptic functions. Therefore the radial ve¬ 
locity behaves as —?> 0 for large r and z/ —>■ 1 for r —>■ 2M, independent of the 
value of C. 

(3) Vanishing Ricci curvature: 

Setting R{h) = 0 in Eq. (fT71) gives r(l — = const, leading to 

1 

V^l -k CN'^rjM ’ 


1/ = 


(29) 
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Fig. 2. Schwarzschild spacetime. The comparative behavior of the slicing function f(r) (left) and 
the corresponding radial velocity u (right) associated with the special observer families considered 
in the text (i.e., CMC, maximal, vanishing Ricci scalar, harmonic and PG) is shown for the 
following choice of parameters: C = 10, k = dil (CMC); C = 10 (maximal); C = 1 (vanishing Ricci 
scalar); C = 4 (harm). The outgoing PG slicing {ly > 0) is shown to compare with the remaining 
mostly outgoing radial coordinate associated observers, apart from the positive curvature CMC 
case which has both radially ingoing and outgoing regimes of their associated observers compared 
to the usual static observers to which all of these are compared (u = 0). Note that / diverges at 
radial infinity for the zero curvature, PG, and constant mean curvature slicings, and at r = 2M 
for all but the constant mean curvature slicings, while the maximal and harmonic slicings are 
chosen to vanish at radial infinity. 


and 

/M = |L(l_„) + 2Ml„(i^) . (30) 

Therefore the radial velocity behaves as —>■ 0 for large r and :/ —)> 1 for 

r —>■ 2M, independent of the value of C. 

(4) Harmonic: 

The harmonic condition (l30l) gives 

0 = drV + — , (31) 

r 

so that V = CM"^/r"^ and 

f{r) = CM\nN. (32) 

C = 4 is a common choice in the literature. Therefore the radial velocity behaves 
as —>■ 0 for large r and 1 for r — 2M (for that choice of the integration 
constant). 

Let us investigate now spherical slicings characterized by intrinsic flatness, i.e., 
by the vanishing of the spatial Riemann tensor. Examining Eq. (US it is sufficient 
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to impose the condition 

L=107=^, u = ±Vl-N\ (33) 

which obviously yields a flat induced metric apparent from Eq. ([8]) since hrr = 
1, while similarly the acceleration vanishes as well from Eq. (HU, leading to a 
geodesically parallel slicing. Because of the above choice of the lapse function, the 
slicing is associated with the “rain” coordinates introduced above, whose adapted 
observers are known as Painleve-Gullstrand (PG) observers [314] . They are geodesic 
and irrotational, and admit orthogonal time hypersurfaces which are intrinsically 
flat. Their radial velocity relative to the static observers is = —\/l — N'^ = 
— ^j2M/r, having selected the minus sign corresponding to a free fall from rest at 
infinity (ingoing PG, for example). Integrating Eq. Q then gives 

/(r) = -2V2Mr + 2M\n ^^ . (34) 

However, PG slicings are not extrinsically flat since 


MK{Af) = 


2M3 


(g) (g) sin^ 0 (g) 


(35) 


PG observers are in a sense “complementary” to the static observers, which instead 
have extrinsically flat orthogonal time hypersurfaces. 


5. Splitting of the curvature tensor by a generic family of 
observers 

For rotating black holes, the situation is more complicated, so we examine first the 
splitting of the curvature tensor using the test observer congruence associated with 
a given stationary spacelike slicing. Let U be the unit timelike 4-velocity vector 
tangent to the test observer world lines orthogonal to the slicing. The orthogonal 
decomposition of the covariant derivative of this irrotational congruence U is 

= -uMuf ( 1 ) 

where a{U)^ = U°‘VaU^ is the acceleration and 

e{uu = P{uzp{uYpV(^u,) ( 2 ) 

is the expansion tensor with trace 0(t/) = 9{U)°‘cn the expansion scalar usually 
referred to as just the expansion of the congruence, while the tracefree part of the 
expansion tensor cr[U)ai 3 describes the shear of the congruence. 

Details of the decomposition of the curvature tensor and corresponding Einstein 
equations may be found in many places, including Ref. [24) . Using an observer- 
adapted frame, the latter lead to the following set of equations in vacuum, setting 
the Ricci tensor to zero 
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0 = R^r = [V{U), + a{U)MUr - V(C/)(ue)0(t/) - Tr [9{Uf], 

o = R^a = -^iu)Mura - e{U)d\] , 

0 = R\ = -[V(t/)(iie) + 0([/)]0(t/)% - [V(C/)(, + a(t/)(b]a(C/)“) + i?(C/)^ ,(3) 

where the index T corresponds to the component tangential to the congruence U 
(and proportional to the orthonormal temporal component). Here V(t7) = P(C/)V 
and V(t/)(iie) = P{U)£u denote the spatial covariant derivative and the spatial 
Lie derivative, respectively. 

Introducing the spatial symmetric “strain” tensor |28i29i:ioi:iii^ . 

SiU)ab = [V(C/)(b + a(17)(b]a(C/),), (4) 

these equations can be written 

V(C/)(Ue)0(C/) = Tr [S{U)] - TT[0{Uf ], 
V(C/)J0(t/)%-0(C/)d%] = O, 

R{urb = [V(t/)(ne) + eiumurb + SiUTb , (5) 

If 17 is a geodesic congruence (i.e., a{U) = 0), the strain tensor then vanishes 
identically and Eq. ([S]) reduces to 

V(17)(iie)0(t/) = -Tr[0(17)2], 

V(17)40(t/)%-0(C/)d“e] =0, 

RiuTb = [V(c/)(ue) + eiUMurt ■ (6) 

Therefore in vacuum for irrotational geodesic slicings one has the general result 
that a nonzero extrinsic curvature implies a nonzero spatial Ricci curvature. This 
is exactly what happens in the Kerr case when considering Painleve-Gullstrand 
observers. What is really special in the Schwarzschild case, instead, is that the 
extrinsic curvature associated with Painleve-Gullstrand slicings is nonzero but 

[V(C/)(ue)+0(C/)]0(G)% = O, (7) 

implying that R{U)°'b = 0. We will discuss these as well as other related properties 
in the next section for a Kerr spacetime. 

6. Kerr spacetime 

Consider the Kerr spacetime representing a rotating black hole with its line element 
written in standard Boyer-Lindquist coordinates (t, r, 0, 4>) as 

ds^ = -(l- dt^ - sin^ 9dtd(j)+^dr'^+J: d9^ + ^ sin^ 9 dcj)'^ , (1) 

where A = — 2Mr + a^, E = cos^ 9 and K= — Aa? sin^ 9. Here 

M and a with |a| < M are the total mass and the specific angular momentum of 
the source respectively. The event horizons are located at r± = M ± y/M'^ — a?. 
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The Boyer-Lindquist t = const slicing of the Kerr spacetime is a well known 
maximal slicing associated with the zero angular momentum observer (ZAMO) 
family of fiducial observers, with 4-velocity field orthogonal to the slicing leaves 
given by 


= -Ndt, n = N-\dt- N^d ^,), (2) 

where N = (—= [AS/A]^^^ and N‘^ = gt,p/g,p 4 , = —2aMr/A are the lapse 
function and the only nonvanishing component of the shift vector field, respectively, 
satisfying the relations (useful in manipulating expressions) 


g<p<pN‘^ = A sin^ 9, 


N‘t> 

JP 


2aMr 

AS 


An orthonormal frame adapted to the ZAMOs is given by 


(3) 


Cf = n , Cf = — - = df , = 


-\/ grr 




— y ^4> ~ 




d4> = di, (4) 


with dual frame 


(jj 


= Ndt, uj^ = y/^rdr, J = ^{dcj> + N^dt) . (5) 


The ZAMOs are accelerated and locally nonrotating, and have a nonzero but 
tracefree expansion tensor completely described in terms of the (shear) expansion 
vector 9{n)^°‘ = 0{n)°‘^ as 

9{n) = e^^e{n)^ + 9{n)^^e^ . (6) 

The extrinsic curvature K(n) = —9{n) is nonzero, whereas its trace vanishes making 
the slicing a maximal one. Tr [Ar(n)^] also vanishes, whereas Tr \K{nY] is nonzero. 
In terms of the dimensionless inverse radius u = M/r and the dimensionless rota¬ 
tional parameter a = a/M, the latter turns out to be 


Tr [K{nY] = 2a^vA sin^ {6) 


. 2— 4/3E-I-8u(l-I-a^M^) 


[SA-k 2u(l-f a2u2)]223 


with 


a = —3 -I- 8u — 6a^u^ + cYu ^, P = —3 -I- 4u — + 2d^u '^, 


2„,2 I oa2„,3 


and 


E = 1 -I- aPP cos^ 9, A = 1 — 2^ -|- . 


(7) 

( 8 ) 

(9) 


In the weak field limit m <C 1 the previous expression becomes 


m 2 Tr [A:(n)2] = 18 sin2 9 


1-i(4 + 13cos2 0)a2u2| +O(ui0). (10) 

o 
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Similarly, the intrinsic curvature of the Boyer-Lindquist time coordinate slices 
is nonzero. In fact, the spatial Ricci tensor has the approximated expression 


(R(h)“b) = 


+ 


-2 0 0\ 

/ 

'0 1 0\ 

0 1 0 -k 90^ 

cos 0 sin flu"* 

0 0 0 

0 

0 

\ 

, 000 / 

/5 cos^(0) — 

1 —6cos0sin0 

0 ' 

0 

1 — 3 cos^ 9 

0 

V 0 

0 

—2 cos^ 9 


0{u^). ( 11 ) 


The three matrices appearing here have zero trace and only at the next order u® do 
nonzero trace terms appear. In fact, R{h) = Tr [K{n)^], which results from 
combining the first and the third of the general relations ([5|) when 0([/) = 0, and 
since Tr [K{n)'^] starts as 0{u^) in its asymptotic expansion (see Eq. (fTOll L 
The ZAMO kinematical quantities only have nonzero components in the r-9 
plane of the tangent space, i.e. 

a{n) = a{ny Cf + a{ny = V(n) In N , 

9{n)^ = 9{n)f e, + 0(n)/‘ eg = -^V(n)iV^ 

Eor later use it is convenient to introduce the quantity 

_ N 

Vn — I — / ! 

as well as the curvature vectors k{x^,n) associated with the diagonal metric coeffi¬ 
cients |24i;i;ii:i4i:i5| 

k{x’',n) = k{x'', nYef + kl^x'', nYeg = — V(n) lny/gY. (14) 

The relevant ZAMO kinematical quantities are listed in Appendix A. A recent 
review of ZAMO slicings can be found in Ref. [36] . 

Let us focus on axisymmetric slices / = f{r,9), because the Kerr metric is 
axially symmetric. The timelike unit normal to these slices is given by 

A/"'’ = —L{dt + frdr + fed9 ), Af = 'y(n + , (15) 

with relative velocity components and associated Lorentz factor 

Va = ^^fa, = {I- . (16) 

y 9aa 

The induced metric on E is given by 

fgrr 0 0 \ / 0 0 /A / /2 Yfe 0\ 

{hab) = 0 gee 0 + 5*0 0 0 fe \ + gtt \ frfe /e 0 , 

V 0 0 g^J \frfe0j Vo 0 0 / 


( 12 ) 

(13) 


( 17 ) 
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or explicitly 


hrr — 9rr ( 1 ^2 ^ ’ ^96 — 996 f 1 ^2 j ’ 


hd,d> — 


hr9 — 


\J 9rr966 


— yj9'^'^9<P<P ; 


IZnVg . (18) 


The new lapse is L = 'jN, since N’^fc = 0, taking into account that Na = N^S'^a, 
whereas the new shift is specified either by the covariant components 


Li — 9ttfr — 


9rr 


In 


t'f , L 2 — 9ttf9 — — 




'§ ! 1^3 - 9t(tl - -^0 J (19) 


or the contravariant components 

Vf 


= -YN- 


= -y^TV- 




( 20 ) 


\/^ \/^ 

Finally one can evaluate the (nonorthogonal) basis Ec and its dual frame 11° using 
Eqs. (HOI), i-e., 

El = dr + frdt, E 2 = dg + fgdt, E 3 = d(ji, (21) 


and 

= dr + ^^Af\ =d0+ = d(^-(22 ) 

V 9rr QQQ y/9^ 

Having written the new spatial metric and the normal congruence Af, obtaining 
both the kinematical quantities of Af as well as the extrinsic and intrinsic curvature 
of the slice E is now straightforward. The 4-acceleration a{Af) = Vfj-Af turns out 
to be 


a{Af) = a{Af)in^ + a{N) 2 ^'^ , (23) 

with 


a{N)i = yf^{dr\y3L) = y/^ [a(n)f +"f'^{vfdrVf + VgdrVg)\ , 

a{Af )2 = yfm{dg InL) = y/^ [a{n)g A- -^'^ivfdgi'f + VgdgVg)] . (24) 

With respect to the spatial frame {Eq}, the components of the extrinsic curva¬ 
ture tensor K{Af) = K{Af)ab^°' ® are given by 


K{Af)ii 
K{Af)i2 
if (AT) 13 


9 rr^ [^f(Wf T ^f|r] ; 

-V9rrg99l I'fiWg - V^Zg) - UgVnZf + Vg 
Va rr 94>4>'yZf , 


K{Af) 

K{Af) 

K{Af) 


22 = 999 '! [l^giWg - 2 VnZg) A Vg^g 

23 = \/g99g<t>4>iZg , 


33 — 


(25) 
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where 


• v) - [9{n)^a “ {9{n)^ ' v)l^a] , 

Wa = n)-v) + [a{n)a - (a(n) • z/)z/a], (26) 


and = P{n)°'aP{nYiV 0 Va, i-e., 

= df-Vr — k{r, n)gVg , 

v§\f. = dfUg + k{r,n)gVf , 

The trace then turns out to be 

Tr [K{N)] = -7 I + -^+1 


— dgi'f + k{0, n^fVg , 
r'eie = dgi^g - k{ 0 , n)fVf , 


^f\f + + I^fWf + I'gWg 




(27) 


(28) 


The intrinsic curvature associated with the induced metric on E can be easily 
calculated too. The nonvanishing components of the spatial Riemann tensor are 
given by 


R{.^f)r4>r(p — grrg<j)(p[—E{n)gg + E{n)^^Vf] + K(Af)uK(Af)33 + [K{Af)i3]'^ , 

R{Af)r4>0<i> = gci,4,^grrgee[E{n)-g + E{n)^^VfVg\ - K{Af)i2K{Af)33 + K{Af)i3K{Af)23 , 
R{Af)rer<i> = grry/geeg<j><j>[{H{n)ff + 2H{n).g - i'nE{n)~g)vf - {H{n)~g - VnE{n)gg)vg] 
-\-K{Ar)iiK{Af)23 + K{Af)i2K{Ar)i3 , 

R(Af)rd9^ — g99'\/ grrg4><p [( E (j^hg T PjiE{Tl)ff')^f P {^H{Ti)gg + 2i7(7r)ff + VjiEijl) ~g^U^ 
—K{Ar)i2K{Af)23 — K{Af)22K{Af)i3 , 

R{Af)9<t>9<i> = g99gti>ti>[—E{n)ff + E{n)^^vf\ + 77(A/”) 22 (A /”)33 + [K{Af)23\^ , 

R{Af)r9re = grrge9{-E{n)^^ + [E{n)gg + Vri{2H{n)f.g - VnE{n)ff)\vl 
-2[(1 + vl)E{n)f.g + Vn{H{n)ff - H{n)gg)\vfVg 

+ [E{n)ff - Vn{2H{n)fg + VnE{n)gg)]vj} 

-K{M)iiK{N)22 + [K{N)i2? , (29) 


where E(n) and H{n) denote the electric and magnetic parts of the spacetime 
Riemann tensor as measured by ZAMOs, respectively, defined by 


E{n)a 0 = Rati/Sun^n '', H{n)a 0 = -[R*]agii3i^n^"n''. (30) 

Their nonvanishing frame components are listed in Appendix A. Note that confor¬ 
mally flat axisymmetric slices do not exist in general, as shown in Refs. Pile] . 

We know of no time slicings which take advantage of the extra freedom to depend 
on the polar angle 0. All of the interesting slicings known to us fall into the special 
case of “spherical” slicings where / = /(r) depends only on the Boyer-Lindquist 
radial coordinate r, so that Vg = 0 and Af = 7(71 -I- Uf i.e., it is just a boost of 
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n in the radial direction. Thus the new observers moving orthogonally to the new 
time slicing follow radial trajectories. The induced metric on E simplifies to 

fgrr 0 0 \ 

(hab) = 0 gee 0 + gt^fr 0 0 0 + gttfr 0 0 0 . (31) 

V 0 0 g^J Vl 0 0/ Vo 0 0/ 

The components of the 4-acceleration of observers having world lines orthogonal 
to the slices E become 


n(A/*)i — ^J~g^{dr In Z/) — ^J~g^ + y , 

a{M )2 = InT) = [a{n)g + , (32) 

while the components of the extrinsic curvature tensor of each slice are 
ZL(A/")ll = grrT [^f(fbr + dr^r\ , 

K{J\f)i 2 = -y/grrgee'll^f a{n)g + i/„6»(n)^g + k{r,n)g , 

-^(•^)l3 “ grrg4>(p^^r ; 

K{M)22 = geeivfk{0,n)f, 

K{Af)23 = y/99e^"f0{n)^g , 

K (A/')33 = g^<j>'^vfk{(j), n)r , (33) 

with 


Tv[K[J\f)] =—^+ {a{n)f — k{9^n)f. — k[(j),n)f.)vf\ . ( 34 ) 

Finally, the nonvanishing components of the spatial Riemann tensor (1291) sim¬ 
plify to 

i?(A/')r0r0 = grrg<t,<t,[-E{n)gg + + K{M)iiK{M) 33 + [KiN')l3]‘^ , 

R{J^)r(j>e<p = g<j><j>\/grrgse E{n)~g — K{M)i 2 K{M )33 + K{M)i3K{M)23 > 

R(/y)r9r4> ~ grry/~g99^(y4>\E(jl)rr + ‘^E(ji)^q b'^E( tI^ 

+K{J\f) II K{J\r) 23 + K(J\f)l2K{J\f)l3 , 

R{J^^r99(fi — g99y/grrgtptp [ Eli^Tl^^g + 

— K{M)i2K{M)23 — K{M)22K{N)i3 , 

R{'N')9<j>94, = —g99g<j><j>E{n)ff + K{M)22K{J\f)33 + [K{J\f)23Y J 

R(J\f)r9r9 — grrgSO^ E{tl)^^ + \E{ti)qq + UjiiOiU (ti) ly^E 

-K{N)iiK{N)22 + [K{M)i2f . ( 35 ) 


Note that again as in the Schwarzschild case, geodesic slicings correspond to 
spatially constant lapse function L = Lq, and one could in principle consider the 
analogous hail, rain and drip coordinate systems. However, in this case the simple 
geometrical properties of the intrinsic curvature are lost, and we will limit our 
discussion (see below) to the case of rain observers {Lq = 1), which are identified 
with PG observers. 
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Fig. 3. Kerr spacetime slicings. The behavior of the linear velocities associated with the special 
observer families considered in the text (i.e., harmonic, PG and RT) is shown for the choice of 
parameters a/M = 0.5 and 9 = 7r/2, for which the outer horizon is located at ^ 1.867M. 


Next some explicit examples are considered. The behavior of the associated 
linear velocities is shown in Fig. |31 


6.1. Harmonic slicings 

The harmonic condition (1301) implies 

dfi^f = [k{0, n)f + k{4>, n)f\vf , (36) 


so that 


or equivalently 


Further integrating then gives 




C sin 9 




C 


— r_ 


In 


r — r+ 
r — r_ 


(37) 


(38) 


(39) 


Regularity of the new lapse function then requires C = + o^. 

Harmonic slicings are neither intrinsically nor extrinsically flat. For example, 
the Ricci scalar has the following approximate expression 


m2 R{h) = ia^u'^lcos^ 9 + (3cos2 9 - l)u] + 0(m®) . (40) 


The traces (linear, quadratic, cubic) of the extrinsic curvature associated with har¬ 
monic observers are shown in Fig. |4l 
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2,4 


2,2 


2,6 


2,8 


r/M 



Fig. 4. Kerr spacetime, harmonic observers. The traces (linear, quadratic and cubic) of the extrinsic 
curvature associated with harmonic observers are shown for the same choice of parameters as in 
Fig. [3] 

6.2. Geodesic slicings 

Geodesic slicings are characterized by vanishing acceleration, i.e., by constant lapse 
function L = Lq oi the spacetime metric written in adapted coordinates. Therefore, 
they are identified by the condition 




(41) 


the same condition as in the corresponding Schwarzschild spacetime. 

When Lq = 1 these slicings are associated with the PG observers, whose world 
lines are a timelike congruence with both vanishing acceleration and vorticity in 
the Kerr spacetime as well |:ii4ii4i:i7i:i8l^ . Their relative velocity with respect to 
the ZAMOs is 



(42) 


corresponding to radially infalling observers, leading to 



(43) 


In contrast to the Schwarzschild case, the geodesic condition does not imply 
intrinsic flatness since the Ricci tensor of the induced metric does not vanish iden¬ 
tically. It can be expressed in terms of the extrinsic curvature tensor as in Eq. 

The spatial curvature scalar is given by 



( 44 ) 























September 17, 2015 0:19 WSPC/INSTRUCTION FILE kextr'final'grqc 


Slicing black hole spacetimes 25 



rIM 


Fig. 5. Kerr spacetime, PG observers. The rescaled traces (linear, quadratic and cubic) of the 
extrinsic curvature associated with the PG observers are shown for the same choice of parameters 
as in Fig. [3| 


Thus the geometry associated with the PG observers in the Kerr spacetime is 
neither intrinsically nor extrinsically flat. 

For completeness we list also the linear, quadratic and cubic invariants of 
K{Nvg)- 


MTr[if(VpG)] 
Tr [KiAfpof] 


Tr [E:(Vpg)^] 


■7/3/2 

-(3 + a^v?). 


S[2(l + a2u2)]i/2 


-[9 + (10 — 3 cos^ 0)a^u‘^ 


2E3(l + a2it2)' 

+(5 — 6 cos^ 0)a‘^u‘^ + cos^ 9 


,9/2 


-[15 + 3(7-008^ 9)a^u‘^ 


E4[2(l + a2rt2)]3/2' 

+21 sin^ 9a^u^ + (7 — 9 cos^ 9)a^u^ + cos^ 9 , 


(45) 


The eigenvalues of iL(A/pG) are also easily evaluated 


MAo = 
MX± = 


V2u^/^ 


tVTT 




MXn 


(1 + a^u^) ± 


(46) 


/ \ 3/2 

j [9 + 3(3 — 5 cos^ 9)a?u^ + cos^ 9. 


Note that Aq is always positive. In the Schwarzschild limit the above expressions 
reduce to MXq = MX+ = —2MA_ = . The traces (linear, quadratic, cubic) 

of the extrinsic curvature associated with PG observers are shown in Fig.jSl 
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6 . 3 . Hyperboloidal slicings 

Numerical relativity computations (like those concerning outgoing gravitational 
radiation) use preferred hyperboloidal spacetime slicings. A general framework for 
the construction of hyperboloidal coordinates with scri-fixing (i.e., with an explicit 
prescription to fix the coordinate location of null infinity) for stationary, weakly 
asymptotically flat spacetimes (including black hole spacetimes) was first developed 
by Moncrief [20] (see also [41142] b 

Recently this method has been successfully applied to the numerical investi¬ 
gation of tail decay rates in a Kerr spacetime by Racz and Toth (RT) [43] (see 
also Refs. [44145146] !. Their construction starts from the standard Boyer-Lindquist 
coordinates {t, r, 9, cj)), passing then to ingoing Kerr coordinates {t, r, 9, such that 

i =t — r + J ^ ^ ° dr , ^ = cj) + J ^dr . (47) 


Finally, the time and radial coordinates t and r are replaced by the new time 
coordinate T and the compactified radial coordinate R, implicitly defined by 


t = T-AM 


In 1 - 


M2 


AM^-R^ 


2M^R 

M^-R^ 


(48) 


The advantage of using RT coordinates (T, R, 9, 0) is that the time slices are horizon 
penetrating and connect to future null infinity, so that no boundary conditions are 
needed. The relation for the radial coordinate r = r{R) can be easily inverted, 
leading to 


R _ r ae m 
M V r2 r 


(49) 


Solving for T then gives T = t — f{r) with 
r^ -I- a 


/(r) = r-J 


■ dr — AM 


In 1 — 


R^ 


'21 


IM^ + R 

A M^-R^ 


(50) 


RT observers thus form an irrotational congruence A/rt of world lines moving ra¬ 
dially with respect to ZAMOs, with a relative velocity given by 


N 

y/ 9rr 


2Mr AM 
A T 


1-f 




R 


M2 


(51) 


with R = R{r) as defined in Eq. (HHl) . It is easy to show also that 


lim Vf = — A , lim izf = 1, (52) 

r—r—>-oo 

with Rf passing monotonically from —1 to -1-1, irrespective of the value of 9 (see 
Fig. [3]). 

The traces (linear, quadratic and cubic) of the extrinsic curvature are shown in 
Fig.d 
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Fig. 6. Kerr spacetime, RT observers. The traces (linear, quadratic and cubic) of the extrinsic 
curvature associated with RT observers are shown for the same choice of parameters as in Fig. |3] 


7. Concluding remarks 

Special congruences of timelike world lines in Kerr spacetime have been studied 
extensively in the literature (see, e.g.. Refs. |33I34] and references therein), while in 
contrast there exist only a few interesting spacetime spacelike slicings that have been 
investigated. The most natural one is the maximal spacelike slicing associated with 
the Boyer-Lindquist time coordinate, whose properties are best described in terms 
of the orthogonal congruence of timelike world lines threading this slicing associated 
with the so-called ZAMO observers. This congruence is accelerated and shearing 
but (locally) nonrotating, but the induced metric on the slices does not have any 
special properties. Nevertheless, other interesting slicings can be found compatible 
with the Killing symmetries of the spacetime, motivated either by geometrical or 
analytical considerations. 

We have developed here a general framework for analyzing all possible stationary 
slicings allowed in a given stationary geometry, then specialized our results to black 
hole spacetimes, reviewing the various familiar examples which can be found in the 
literature. Geometrical properties of these slicings can be characterized either in 
terms of those of the kinematical properties of their orthogonal timelike congruences 
interpreted as the world lines of a family of stationary observers, or in terms of 
intrinsic or extrinsic curvature properties of the slicing itself. We have considered 
in detail the case of slicings associated with a new time coordinate function of the 
Boyer-Lindquist coordinates {t,r). The case of geodesic observers includes the so- 
called Painleve-Gullstrand slicing of Kerr. The latter are even more special in the 
limiting Schwarzschild case where the slices are also intrinsically flat. Finally we 
have also presented examples of “analytical” slicings of black hole spacetimes, like 
those associated with harmonic and hyperboloidal time gauges. 
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Appendix A. ZAMO relevant quantities in the Kerr spacetime 


We list below the non-vanishing components of the electric and magnetic parts of 
the Riemann tensor as well as the relevant kinematical quantities as measured by 
ZAMOs. 

The components of the acceleration and expansion vectors are given by 

ainf = {-2r\r^ + a^)A + _ 4Mr^]} , 

, 2a^Mr ,9 9 , . ^ 

a{ny = -^ 7 --(r^ -I- a^) sim 9cos9 , 

94,4> 

/ , lUnSm ^r^9/9 9\ 1/9 9\i 

[ 2 ^ + a ) + - a )] > 




2a^Mr 


Va cos 9 sin^ 9 , 


and the curvature components are 
E(r — M) — rA 


k{r,ny = 


E3/2yA 

r-\/A 


k{r, n)^ = 


sin 9 cos 9 


S3/2 


sin 9 cos 9 


k{9,ny = -^^, k{9,n) = -^ 


k{(t>,ny = 


•\/Asin 9^ i\/r\ A/rx'(oJ 2 , , „ 2 \ 




■[-E2(r - M) - ME(3r2 -h a^) + 2My{y + , 


k{(t),ny = - 


cos 9 sin 9 


[E^A + 2 Mr(r^ -f a^]. 


(All) 


(A.2) 


Finally, the nontrivial components of the electric and magnetic parts of the Riemann 
tensor with respect to ZAMOs are given by 

Eff. = dfa{n)f + a{n)f — 30(n)?^ — a{n)gk{r, n)^ , 

Egg = -Eff + dfk{9, ny - [k{9, n)’’]^ -f dgk{r, ny - [k{r, n)^]^ , 

E~g = dga{n)f + a{n)g[a{n)f + k{9, n)’'] - i0{n)f^9{n)g^ , 

Hff = dg9{n)~^ - 9{n)f,^k{(l), nf + 9{n)g^[k{9, n)’’ - ny ], 

H-g = -df9{n)~^ + 29{n)f.^k{(t>, n)’’ -f 9{n)g^k{r, nf , 

= ^(n)~^[a{n)g + k{<p, nf] - 9{n)g^[a{n)f + k{<l), n)^] , (A.3) 

in terms of the above kinematical quantities. 

In the equatorial plane, the acceleration and expansion and curvature vectors 
are all directed along the radial direction, with components 

i \r _ ^ + yy — Aa?Mr n ( \r _ aM(3r^+a^) 

“ r^^A r3 + a2r + 2a2M ’ ~ ~ + a^r + 2a‘^M) ’ ^ ’ 
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and 


n)’’ 

k{4>,ny 


Mr — 
r'^yfK 


k{9, n)’’ 


(r^ — a'^M)'/A 


^ 2 (j ,3 _|_ g 2 j. _|_ 2a?M) 


J.2 ’ 


(A.5) 


respectively, whereas the nonvanishing frame components (IA.3I) of the electric and 
magnetic parts of the Riemann tensor reduce to 


M(2r^ + 5r^a^ — 2a^ Mr + 3a^) 

j,4(j,3 _|_ q 2 j , _|_ 2a^M) 

SMa^"^ + a'^)\/A 

^4(j,3 _|_ q2^ _|_ 2a^M) 


^ee ~ ^rr , 


~ 


M 

^3 ’ 


(A. 6 ) 
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